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ABSTRACT 


This report covers several areas that were studied during the 
year, all of which concerned observations that were made or could be made 
with the radio telescope at the La Posta Astrogeophysical Observatory 
of the Naval Electronics Laboratory Center in California. In the area 
of resolution enhancement, the use of Fast Fourier Transform programs 
was investigated for possible application to millimeter wavelength maps 
of the Sun. A difficulty that arises with the La Posta maps is that 
they are limited to 35 arc-minutes square while the smeared out solar 
image is larger than that. 

A list of possible cometary emission lines near 13 millimeters 
is presented. Although preparation of the list was inspired by the 
appearance of Comet Kohoutek, the results are applicable to any future 
comet. 

The brightness temperature of the Sun at 8.6 millimeters was 
measured using the Moon as a calibration source. The result does not 
confirm a deep absorption feature as apparently observed by earlier 
workers. 



N. E. Chatterton, Ph.D. 
Manager, Research Laboratories 
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1. INTRODUCTION 


The acquisition of large quantities of observational data on the 
Sun from the Apollo Telescope Mount on the Sky lab spacecraft and from 
earlier satellites such as Orbiting Solar Observatory 7 has reemphasized 
the importance of using combined approaches and methods of observation 
to study physical phenomena. In other words* data obtained from space 
should be accompanied by ground-based observations in order to accomplish 
the best and most broadly based analyses. Because the ATM and 0S0 experi- 
ments involve primarily the ultraviolet and X-ray regions of the spectrum* 
they are well complemented by ground-based radio observations. The Space 
Sciences Laboratory at the Marshall Space Flight Center has performed or 
has had access to such observations in the millimeter wavelength region, 
and it is important to use these data in studies of the Sun. 

Several areas involving radio astronomy of the Sun and other bodies 
have been pursued and are discussed in this report. These include reso- 
lution enhancement of maps of the Sun, possibly detectable line emissions 
in the millimeter spectral region from comets, and the observation of the 
brightness temperature of the center of the Sun at 8,6 millimeters and its 
comparison with models. 



2. RESOLUTION ENHANCEMENT 


An image can usually be considered to be a convolution of the 
object with an antenna power pattern or point spread function plus noise. 
Normally, filtering of noise and deconvolution of the observed image 
distribution to obtain the true object distribution are accomplished by 
using Fourier transform techniques. Because the millimeter-wavelength 
radio maps of the Sun are usually in the form of digitized arrays, we 
use the finite discrete Fourier transform, whose properties are conveniently 
summarized by Gentleman and Sande (1966) . [We actually use the normaliza- 
tion and some of the conventions given by Brault and White (1971).] 

The discrete Fourier transform is given by 
N 

F(t k ) .= N" 1 2 ^ F ( x j) ex P ( -i 2 tt x j t k ) (1) 

j=l 

and the inverse by 
N 

F ( x j) 5 ^ F (t k ) ex P 2tt x j ^ (2) 

k=l 

where F(x^) are the data, - (j - l)Ax are the points at which the 
data are sampled, Ax is the sampling interval, N » the number of data 
points, and t^. = (k - l)/(NAx) are the discrete spatial frequencies. Use 
of the Fast Fourier Transform algorithm, described, for example, by 
Gentleman and Sande (1966) and by Cochran et al (1967) , has significantly 
reduced the time needed to compute Fourier transforms, thereby making 
them convenient to use. 

The power spectrum is easily computed from the Fourier transform 
by the relation: 

p (t k ) = F(t k ) F(t k )* (3) 

where F(t k )* is the complex conjugate of F(t k ). The power spectrum is a 
convenient quantity to use in examining the influence of noise on the 
data. It will be formally used in the optimum filter. 
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The distribution of intensity or radiance in an image I(x) can 
usually be written as 


I(x) = S(x) 4- N(x) 

W 

00 

= f 0(x f ) A(x - x’) dx’ + N(x) 

(5) 


— 00 


or in the discrete form 

I(Xj) = S(xj) + N(xj) 

N 

= ^ tKxi.) A(xj 
1=1 

where 0(x) is the true object distribution, A(x) is the antenna power 
pattern or point spread function, N(x) is the noise, and S(x) is the 
smeared signal which would be observed in the absence of noise. (The 
two uses of the symbol N should cause no confusion.) The discrete form 
of the convolution shown in Equation 7 must be considered cyclic, i.e., 
the N data points repeat themselves. The cyclic property introduced by 
the sampling process and the finite number of points require care in 
the restoration process. The two ends of a set of data must smoothly 
approach the same value. 

The attempt to recover the true distribution 0(x) form the 
observed distribution I(x) proceeds as follows. The Fourier transforms 
of Equations 6 and 7 are 

Ut k ) = S(t k ) + N(t k ) (8) 

and 

T(t k ) =■ 0(t k ) • A(t k ) + N(t k ) (9) 


x*) + N(x ) 


( 6 ) 


(7) 
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where we use the property that the transform of a convolution is equal 
to the product of the transforms. If we were able to observe S(x), then 
the solution for 0(x) would simply be the inverse transform of 0(t) given 
by 

0(t) = s'(t)/A(t) . (10) 

This process is limited only by the fact that A(t) must normally not 
be zero. If we let 


Ro(t) 


{s’(t)/A(t)'- , A(t) * 0 
(0 ,“A(t) = 0 


( 11 ) 


then Rq(x) is known as the principal solution (Bracewell and Roberts, 
1954). Normally A(t) will be zero beyond some frequency corresponding 
to the Rayleigh limit. In the presence of noise, the simple application 
of the restoration process to the image yields 


I(t)/A(t) 

(12) 

S(t)/A(t) + N(t)/A(t) 

(13) 

V(t) + N(t)/A(t) . 

(14) 


(Again, we have to be careful about 'K being zero.) Because A(t) becomes 
small at high frequencies, Equation 14 amplifies the high-frequency noise 
with disastrous effects on the restoration. To avoid the amplification 
of the noise, an optimum filter is applied (Brault and White, 1971): 


R 2 <t) 


<"*w 

Kt) 

A(t) 


vt) 

p s < t > + P N(t) 


(15) 


where P g (t) and Pjj(t) are the power spectra of the smeared signal and the 
noise. These power spectra can often be modeled after examination of the 
power spectrum of the image. Usually P N can be considered to be a con- 
stant equal to the mean of its high-frequency component. In other words, 
all of the high-frequency components of the observed image are noise, 
and the noise power spectrum is independent of frequency. 
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The development described above will not allow recovery of the 
object distribution at frequencies above that corresponding to the minimum 
separation given by the Rayleigh criterion. In recent years, there has 
been considerable effort devoted to attempts at resolution beyond the 
Rayleigh limit. Several schemes have been devised to extrapolate the 
spectra to higher frequencies. These schemes are based on various 
properties, such as the fact that the true object distribution must be 
positive or that it has a finite extent. However, it has been shown, 
e.g., by Rushforth and Harris (1968) and Miller (1972) that the methods 
based on the finite extent of the object are severely limited by noise. 

The signal- to-noise ratio required to improve even slightly on the 
Rayleigh limit is extremely high; thus, the Rayleigh criterion represents 
a practical limit to what can be realized at reasonable signal-to-noise 
ratios. 

This discussion has been given in terms of one-dimensional 
distributions. Extension of two dimensions is usually straightforward. 

The application of restoration techniques to the 8.6-millimeter 
maps of the Sun obtained at La Posta is subject to several limitations. 

The presence of noise will cause difficulties. One of the main sources 
of noise is the variation of the transmission of the Earth’s atmosphere 
during the time required to make a complete map (about one hour). Another 
limitation is due to the truncated size of the maps. They normally 
cover a solid angle about 35 arc-minutes square, whereas the Sun is 32 
arc-minutes in diameter and the antenna pattern has a half-power width of 
several arc-minutes. Use of the Fast Fourier Transform on the sharply 
cutoff image would create false high-frequency components in the transform 
(aliasing). A possible way around this problem would be to artifically 
extrapolate the image according to a known mean distribution. Of 
course, no details near the limb of the Sun would then be recoverable. 
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3. LINE EMISSION FROM COMETS 


The discovery of Comet Kohoutek more than nine months before its 
perihelion passage in late December 1973 made it possible for observers 
to have considerable time to prepare special equipment and to plan observ- 
ing programs. At that time it appeared that support might be available ' 
to build a line receiver to be used on the 60-foot antenna at the La Posta 
Observatory of the Naval Electronics Laboratory Center near San Diego, 
California. Because it would be desirable to search for line emission 
from comets, especially from certain molecules, an effort was made to 
identify which lines might be detected. Before listing the possible 
lines, we present a brief summary of the characteristics of comets, based 
in part on the review by Brandt (1968) and the proceedings of the Tucson 
Comet Conference (Kuiper and Roemer, 1972). 


3.1 PHYSICAL CHARACTERISTICS OF TYPICAL COMETS 

Comets are usually described as having three parts: a solid 

nucleus, a coma of dust and neutral molecules, and a tail of ionized 
molecules or dust. 

The nucleus is thought to be 1 to 10 kilometers in radius and 
composed of an "icy conglomerate" (Whipple, 1963) containing frozen com- 
pounds such as H^O, NH^, CH^, CO^, etc * Embedded in the icy material 

is meteoric material. The mass may be on the order of 10 16 to 10 2 1 grams. 
The nucleus has never been directly observed. 

The ices in the nucleus presumably sublimate in vacuum at temper- 
atures on the order of several hundred degrees Kelvin and dissociate to 
form the simpler daughter molecules which are observed in the coma. The 
parent or precursor molecules themselves, such as H^O, NH^, and CH^,, have 
not been observed but are inferred from the existence of the observed 
daughter molecules such as the radicals OH, CH, NH, and NH^. Because the 


6 



parent molecules have not been observed, it is, therefore, important to 
search for them in order to better understand the chemistry of comets. 

The radius of the coma is typically about 10 3 kilometers, corresponding 
to more than 2 arc-minutes at a distance of one astronomical unit. There- 
fore, the coma would make a suitable target for observation in the milli- 
meter range by the La Posta radio telescope. However, the parent molecules 
would probably exist only in a much smaller volume of the coma surrounding 
the nucleus. Typical expansion velocities of molecules ejected from the 
nucleus are on the order of 1 km/sec. Because the total gas density is 
low, ranging from about 10 12 to 10 14 molecules /cm 3 near the nucleus to 
IQ 2 to 10 4 molecules/cm 3 at the outer boundary, association processes are 
not important. In other words, once the parent molecules are dissociated, 
the daughter components remain apart and never recombine to form the 
parent. Dust and metallic atoms are also observed in the coma. A fairly 
recent development has been the observation of a hydrogen Lyman— a halo 
which extends out to a radius of about 10 6 kilometers. 

The third major element of a comet is the tail, reaching up to 
10 B kilometers or one astronomical unit in length. There are two types 
of tails; occasionally both types appear in the same comet. Type I tails 
are composed of ionized molecules and atoms with some neutral atoms and 
are straight, pointing, to a first approximation, away from the Sun. 

Type II tails are strongly curved and are composed of dust showing a 
reflected solar spectrum. 

In the visible spectrum and for hydrogen Lyman-a, the predominant 
excitation mechanism is resonance fluorescence by sunlight. In the case 
of the forbidden oxygen line, the atom is presumably formed in the excited 
state during the photo-dissociation of its parent molecule. It is not 
known what the excitation mechanism for microwave lines would be — 
radiative, eollisional, or excitation — during the formation or release 
of the molecule. Relative populations of levels emitting the observed 
lines could show how the levels are excited. If they are collisionally 
excited, densities or limits on densities may be inferred. 
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3.2 


OBSERVATIONAL POSSIBILITIES 


The spectral region around 13 millimeters appeared to have several 
lines of interest, and as there would not have been any technical difficulty 
in building a receiver for that wavelength, only lines in that spectral 
region were considered. Table 1 contains a list of lines near 13 
millimeters which might be observed in comets. Most of the information 
on molecular lines is based on the Microwave Spectral Tables (National 
Bureau of Standards Monograph No. 70). After the information was compiled, 
it was discovered that Huebner (1970) had made a more complete listing. 

Also shown in the table are some hydrogen and helium lines taken 
from Lilley and Palmer (1968) . These lines are known as recombination 
lines because they normally are emitted following recombination of an ion 
with an electron. Because the densities are so low* recombination is - 
likely to be rare and the lines are probably too weak to be observed. 

The coma would provide the best location to search for line 
emissions. 

Unfortunately, no financial support became available to enable 
observations to be made of Comet Kohoutek. However, it will still be 
desirable to search for line emissions in the microwave spectra of future 
comets . 
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TABLE 1. MICROWAVE LINES WHICH MIGHT BE OBSERVED IN COMETS 


MOLECULE 


• 

FREQUENCY 


OR ATOM 

TYPE OF LINE 

TRANSITION 

(MHz) 

REMARKS 

H 2 O 

(Water) 

Rotational 

5 23 " 6 1 $ 

22235.22 

Levels are 446 cm - 1 
above ground state 
(equivalent to 22.4 pm). 

nh 3 

Inversion 

( 1 . 1 ) 

23694.49 

( 6 , 6 ) level is almost 

(Ammonia) 

Doubl et 

( 2 , 1 ) 

23098.78 

300 cm - 1 above ground 



( 2 , 2 ) 

23722.63 

state; others are all 



(3, 2) 

22834.17 

lower. All but last two 



(3, 3) 

23870.13 

transitions have been 



(4, 4) 

24139.41 

observed in interstellar 



( 6 , 6) 

25056.02 

medium. 



(5, 5) 

24532.98 




( 4 , 3 ) 

22688.29 


OH 

Lambda 

2 n Vz 9/2, 4 


Levels are about 

(Hydroxyl ) 

Doubling 


375 cm -1 above ground 



F 1 = 4 F" = 5 

23806.5 

state. 



4 4 

23818.18 




5 5 

23826.90 




5 4 

23837.8 


NH 


F 2 (J = 1)-F X (J = 2) 

25586.4 

Predicted by Kerns and 
Duncan 1972. 

HNCQ 


Oo c ~ 1 0 1 

21982. 


(Isocyam'c Acid) 





CH 3 OH 


(J, K) = 4, 1-4, 2 : 

24933. 


(Methylalcohol ) 


5, 1-5, 2 

24959. 



6 , 1 - 6 , 2 

25018. 




7, 1-7, 2 

25125. 




8 , 1 - 8 , 2 

25294. 


H 

Recombination 

6 6 a 

22364.17 


(Hydrogen) 


65a 

23404.28 




64a 

24509.91 




63a . 

25686.28 


He 

Re combi nation 

66a 

22373.28 


(He! ium) 


65a 

23413.82 



64a 

24519.89 




63 a 

25696.75 





4. SOLAR BRIGHTNESS TEMPERATURE AT 8.6 MILLIMETERS 

Several reviews of the solar millimeter spectrum (e.g., Linsky, 
1973; Kislyakov, 1970; Shimabukuro and Stacey, 1968) have indicated 
the presence of an absorption feature or depression in the solar spectrum 
in the vicinity of 6 to 9 millimeters. Because such a feature is not 
expected according to our present understanding of the solar atmosphere, 
it is important to verify its existence. We therefore decided to remeasure 
the brightness temperature of the Sun by carefully determining the ratio 
of the brightness of the Sun and Moon and using the Moon as a standard 
source. This section of the report reviews the background for the 
reasoning behind the measurement, describes the computation of the 
topocentric coordinates of the Moon which were needed, and gives our 
results in which we did not confirm the presence of a deep absorption 
feature. A report on our work, which was carried out jointly with Max 
Bleiweiss of the Naval Electronics Laboratory Center and Fred Wefer of 
the Megatek Corporation, was presented at the meeting of the Solar 
Physics Division of the American Astronomical Society held in Honolulu, 
Hawaii, in January 1974. The abstract of the paper is reproduced in 
Appendix A. 

4.1 BACKGROUND 

Most models of the solar chromosphere, whether of one or more 
components, have the temperature increasing outward. Therefore, as one 
observes at longer and longer wavelengths in the millimeter region, 
corresponding to higher opacities, one should see a monotonically increas- 
ing brightness temperature. However, many reviews of the solar millimeter 
spectrum show a depression or absorption feature. For example, Shimabukuro 
and Stacey (1968) and Kislyakov (1970) plot many observations which have 
considerable scatter but appear to have a minimum near 6 millimeters. 

The exact position (and possibly even the existence) of the minimum 
is poorly defined because of the scatter. On the other hand, Linsky (1973) 
has accepted only those measurements made relative to the Moon and has 
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recalibrated them using his own review of the lunar spectrum. His plot, 
which is shown in Figure 1, shows a depression of over 2000°K at the 
two data points at 8 and 8.6 millimeters. The observation at 8^6 millime- 
ters by Ulich, Cogdell, and Davis (1973) is claimed by them to have 
been verified in unpublished work by Ulich. Published explanations that 
have been offered for the minimum include a multicomponent atmosphere 
with the relative volume of hot material decreasing with height (Zheleznyakov, 
1964) and scattering of the radio waves by plasma- turbulence waves (Kaplan 
and Tsytovich, 1967). 

In view of the importance of irregularities in the millimeter 
spectrum to understanding the details of chromospheric structure, it was 
decided that the existence of the minimum should be verified by a 
measurement at La Posta of the brightness temperature at 8.6 millimeters. 

4.2 TOPOCENTRIC COORDINATES OF THE MOON 

Because the Moon would be used as a standard calibration source, 
it was necessary to be able to point the La Posta radio telescope accurately 
at the center of the Moon. Computer programs already existed at La Posta 
to use astronomical positions as tabulated in the American Ephemeris and 
Nautical Almanac . The tabulated positions are relative to the center 
of the Earth. The Moon, however, is so close to the Earth that parallax 
affects its apparent position as seen by an observer at an arbitrary 
point on the surface of the Earth. Therefore, it was necessary to write 
a computer program to convert from geocentric to topocentric coordinates. 

The relationships between the geocentric and topocentric coordinates 
are summarized below with the particular form and notation taken from 
Raine (1970). The topocentric hour angle h T and declination are given 
by 


sin hrp - A/S 

(16) 

cos h-p - B/S 

(17) 
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SOLAR BRIGHTNESS TEMPERATURE 
{UNITS 1000°K) 


FREQUENCY (GHz) 


300 200 100 80 60 40 30 20 10 



FIGURE 1. OBSERVED SOLAR MILLIMETER SPECTRUM. Points and 
fitted curves taken from Linsky (1973) except 
for A from this investigation. 



and 


sin 6 t = C/F (18) 

where 

A - cos 5 q sin hg (19) 

B = cos <5 q cos hg - p cos <p ' sin n G (20) 

C m sin 6 q - p sin <{> f sin (21) 

D = A 2 + B 2 (22) 

F = (D + C 2 )^ 2 (23) 

and the geocentric hour angle hg is given by 

ti G = ST 0 + 1.00273791 t - A - a ” (24) 


The geocentric right ascension a^, geocentric declination Sq, and 
horizontal parallax of the Moon JI G are tabulated; also tabulated is STq , 
the Greenwich sidereal time at 0 hours Universal Time (U.T.). X is the 
longitude, p is the distance from the center of the Earth in equatorial 
Earth radii, and 1 is the geocentric latitude of the observer, t is the 
Universal Time and 1.00273791 t is the sidereal time elapsed since 0 hours 
U.T. 

Also convenient for the La Posta observers are the hour angles 
from the Greenwich meridian: 

h T,Grw = h T + X ( 25 ) 

and the topocentric right ascension: 

a T = ST 0 + 1.00273791 t - X - h T . (26) 

The La Posta position is specified by 

X = 116°4512 
p - 0.99902436 
<$>' = 32° 30 f 77925 . 
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A copy of the computer program is contained in Appendix B. 

Because the observations were to be made in October 1973, certain parts 
of the program were written with that knowledge incorporated. Use of the 
program for a different time would require slight modifications. The 
input data needed, which are all found in the American Ephemeris and 
Nautical Almanac , are the Greenwich sidereal time at 0 hours U.T. for 
each day and the geocentric right ascension, geocentric declination, and 
horizontal parallax of the Moon for each hour. 

4.3 OBSERVATIONAL RESULTS 

The observations of the Sun and Moon were obtained with the 
60-foot (18.3-meter) diameter radio telescope at the La Posta Astrogeo- 
physical Observatory of the Naval Electronics Laboratory Center near 
San Diego on 26 October 1973. Preliminary measurements were made on 
24 October but were not used in the final results. The date was chosen 
to be near new moon so that the Sun and Moon would be close together in 
the sky, minimizing differential atmospheric effects. The Moon was used 
as a standard source, following Linsky's (1973) suggestion. Measurements 
were made at 8.6 millimeters (35 GHz), at which wavelength the antenna 
has a half-power beamwidth of 2.8 arc-minutes. There was no solar 
activity present near the center of the Sun. 

Two types of observations were made: analog drift scans across 

an entire solar or lunar diameter and digitized small maps of the central 
regions of the disks. The maps normally consisted of 5 x 5 arrays 
with a grid spacing of 1 arc-minute. Table 2 lists some information 
about the scans and maps. Because of the ease of their use in quantita- 
tive work, only maps, particularly one map each of the Sun and Moon, 
were used in the final results. 

The observational process included a measurement of the brightness 
temperature of the sky nearby and at the same elevation as the Sun or Moon. 
The sky temperature was subtracted off in the reduction of the data. 
Internal calibration was carried out by the use of a 3000°K noise tube 
and a 375 °K oven. 
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TABLE 2. TIMES AND POINTING INFORMATION ON DRIFT 
SCANS AND MAPS 
26 OCTOBER 1973 


SOURCE 

TIME 

(Universal 

Time) 

AZIMUTH 

(deg) 

ELEVATION 

(deg) 

Drift Scans 




Sun 

1835 

161.1 

42.7 

Moon 

1912 

167.3 

36.8 

Moon 

2200 

213.2 

30.0 

Sun 

2221 

229.5 

29.8 

Maps 



.. . 4 . . ... 

Sun* 

'1710 


'30.0 

Moon 

1924 

170.7 

37.2 

Sun 

2130 

217.6 

36.2 

Moon 

2142 

209.1 

31.8 

Sun 

2232 

232.0 

26.8 


*Taken from center of large map on the entire Sun. 


Because of the ease of their use in quantitative work, only the 
reduced observations from the maps are presented here. The results in 
units of antenna temperature are shown in Table 3- The source temperature 
is the mean of the 25 or 21 individual data points with the sky already 
removed. The rms deviation a is effectively equal to the standard 
deviation of a single measurement from the mean; it has not been reduced 
by the square root of the number of points. 

The two maps which were ultimately used were selected because they 
were taken most closely together in time — the solar map at 2130 U.T. 
and the lunar map at 2142 U.T. A flaw in the computerized data recording 
process limited the solar map to 21 points instead of 25. The Sun-to-Moon 
ratio is then 
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TABLE 3. ANTENNA TEMPERATURES FROM MAPS 
26 OCTOBER 1973 


SOURCE 

TIME 

(Universal 

Time) 

NUMBER OF 
DATA 
POINTS 

ANTENNA TEMPERATURE 

(°K) 

SKY 

SOURCE 

(Sky Removed) 

a = rms 
DEV. 

Sun 

-1710 

25 


4296 

23.1 

Moon 

1924 

25 

15.9 

102.2 

1.75 

Sun 

2130 

21 

14.1 

4376 

25.8 

Moon 

2142 

25 

16.7 

104.9 

1.76 

Sun 

2232 

25 

26.0 

4335 

19.1 


Is = 4376 

T m 104.9 


41.72 


(27) 


If the brightness temperature of the new moon is taken to be 219 °K 
(Linsky, 1973), then the ratio of brightness temperature to antenna tempera- 
ture is 



219 

104.9 


2.088 


(28) 


The solar brightness temperature is then 9140°K, and the measured rms 
brightness temperature deviations are o s = 53.9°K and o M = 3-67°K. To 
arrive at an uncertainty combining both the solar and lunar measurement 
uncertainty, we use the expression 



to obtain a = 162°K. 


(29) 


The final result for the brightness temperature of the center of 
the solar disk is then 


T b = 9140 ± 160°K 
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The uncertainty is a statistical measure equal to the combination of the 
rms deviations of both the solar and lunar measurements. Of course, 
there are additional sources of error arising from the effects of the 
Earth’s atmosphere, weak activity on the Sun, and irregularities in the 
surface brightness near the center of the Moon. 

Our result is shown as the triangular shaped point in Figure 1 
where it may be compared with Linsky’s (1973) fitted curves — 9641°K at 
8.6 millimeters from the parabolic curve — and with the observed point 
shown as a cross by Ulich, et al (1973) — 7410°K as recalibrated by 
Linsky. Although the point is still below the curves that Linsky fitted 
to the data, it clearly does not confirm a deep depression. Our measure- 
ment is also shown in Figure 2 where, along with Linsky* s curve, it is 
compared to predictions of two recent models — the Harvard— Smithsonian 
Reference Atmosphere (Gingerich, et al (1971)] and a model by Vernazza and 
Noyes (1972). It is notable that, especially at longer wavelengths, 
Linsky’s curve is significantly higher than the model predictions. 
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APPENDIX A. ABSTRACT OF PRESENTATION 


< 


The abstract of presentation at the meeting of the Solar Physics 
Division of the American Astronomical Society held in Honolulu, Hawaii, in 
January 1974 is presented below. The abstract was published in the Bulle- 
tin of the American Astronomical Society, Vol. 6, No. 2, Part II, p. 288. 

Brightness Temperature of the Sun at 8. 6 
Millimeters - Presence of an Absorption Feature? 

W. HENZE, Teledyne Brown Engineering , M. BLEI- 
WEISS, Naval Electronics Lab. Center , & F. WEFER, 

Megatek Corp. - Reviews of the solar spectrum at 
millimeter wavelengths (J. L. Linsky, 1 973, Solar Phys. 

28, 409; A. G. Kislyakov, 1971, Soviet Phys. Uspekhi 13_ t 
495) together with recent work by B.T . Ulich ,(U. Texas,... 
unpublished) indicate the presence of a depression or 
absorption feature near 8 mm wavelength. Because such 
a feature is not expected according to our present under- 
standing of the solar atmosphere, it is important to 
verify the existence of the feature. Therefore, we have 
observed the brightness temperature of the center of the 
solar disk relative to the center of the new moon at 8. 6 
mm with the 60 ft antenna of the Naval Electronics 
Laboratory Center at La Posta, California. Our pre- 
liminary results based on the temperature of the new 
moon given by Linsky (op. cit. ) do not confirm the 
presence of a deep absorption feature or depression 
near 8 mm. 
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APPENDIX B. COMPUTER PROGRAM USED TO COMPUTE 
TOPOCENTRIC COORDINATES OF THE MOON 
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TPCORD - EFN SOURCE STATEMENT - IF N C S > - 


C PROGRAM TPCQRD 

C W. HEN^E OCTOBER 1973 TSE 

C COMPUTES TOPOC.ENTRIC COORDINATES OF LUNAR CENTER 

DIMENSION lbYSID(15) f SDTGT(15)» IDYMN (126),IHRUT(126) ,HAT ( 126 ) 

-1-tHAT Cj 12 6 ) , DC T ( 1 26 ) # RAT I 12 6), 

RHSNPP = .5368C777 "" “ 

RHCSPP , = .842547975 _ 

PI = 3 . 1 A 1 5926536 

XLON.G. * 2.032456858 

C...XLCNG = LONGITUDE IN RADIANS 

WRITE {.0,212 ) 

212 FORMAT ( 1H1, 3 9H SIDEREAL TIME AT GREENWICH AT 0 HR U.T.//1H , 

1 7HCC TUBER/ 1H ,5H OA Y * 5Jh 27HH0UR MIN SECOND RADIANS \ 

DO 1 1 = 1,15 

RE AD( 5 , 111) IDY$ ID< I ) , IHST,_M$T, SST 

111 FORMAT (13, 14, 13 , F7. 3 ) 

-HSJL= IHST 

XMST = MST 

SDTC = HST + UMST + SST/60* )/60. 

SDTGTI I ) = SDT (PI/12.) 

C..,SDTGT = SIDEREAL TIME . .AT GREENWICH AT..0 HR_U,T*_ IN RADIANS 

WRITE! 6, 211) IDYS ID( I ) , IH ST , MS T , SS T , SDTGT U ) 

2.11 FORMAT.! 1H 1 5 ,18, l iu.FJ? , 3 , FI 0 ._51 

1 CONTINUE “ 

WRITE (6, 2 131 

213 FORMAT { 1H 1 , 109H ST MN UT PD T GEQCENT RIGHT ASCNSN GEOCENT DECLINATION 

1ATI0N KQ.R2N.TL PARALLAX ._. TOP DEC..IOPO HA TQPO..HA TGPG..RA/1H , 

2 107H DY DY HR HR HR Mft SEC RAD DG AM ARCSC RAD AM 

3ARCSC RAL RAC LOG ,RM) GRW,RA£ UADI 

DO 2 J=l, 126 

RE AD ( 5 , 1 12 1 10 YMN { J ) , I HRl T ( J ) , I RAHG » IRAMG, R ASG , I CCDG , IDCMG , OCSG 
1 1 1 PRMG t PRSG 

112 FORMAT (13,13, I A, I 3 , F 7 * 3,1 4, I 3 , F 6 • 2 , 1 4 , F 6 • 2 ) 

HRUT = IHRUT(J) 

RAHG. = I RAH G _ 



t ) 


RAMG = IRAMG 
DCOG =: ILCUG 
CCNG = IDCMG 
P.RMG = iPRfcG 

RA = RAHG+(kAHG+RASG/6Q. )/60. 

RAG..jt . kM i P XL 1 2.1 

..•RAG - GEOCENTRIC RIGHT A SC £N5 f UN C F M OON IN K & 0 1 A N $ 
I DSN = 1 ' 

1FU-30) 7,11,11 
11 I DSN = -1 
CCCC = -DCOG 

7 DEC =. UCDC+ ( DCHG+OC SG/ 60 , ) / 6 0. 

IF(IDSN) 8 » 9 t 9 i 

8 DEC ■= -DEC 

9 OCG = DEC*! PI/180. ) 

C .... DC G G.GU.C.E N JR 1C DECLINMI.CN OF MG UN IN RADIANS 
PR = PRMC+PRSC/60. 

- PRG... = PR* CP 1 / LCaCO. I ' 

tv C...PRG= HORIZONTAL PARALLAX IN RADIANS 

N .. TR = 1 • 00 27 3 79 1 * BRUT *(PI/12.) 


TPCORD - .FN SOURCE STATE ME N T - IFN(S) - 

ID - J+8 
ID = ID/9 

SCO = SUTGTI IDJ+TR 

C..-SUG - PRESENT GREENWICH SI CEREAL TIME 
HAG = SDG-X LDNG-R AG 
SNPRG = SIN1PRG) 

CSDCG ^ CCS (OCG) 

SNDCG = S I N { DC G I 
SNHAG ^ SIN (HAG I 
CSHAG = COS ( HAG ) 

A = CSDCG*S NHAG 
B = CS DCG*C$HAG-SNPRG*RHC SPP 
C * SNpC G-SNPK G*RHSNPP 
0 = A * A + B * D 



S = SQK1 (D) 

F = SGRT(UfC*C> 

SNDCT = C/F 

DC rt J) _ARSIN( SNDCT) IlM 7094 

C * * . DC T = T GPOCLN TRIG DEC L I NAT 1 ON IN R AD I AMS 
SiNHT - A/ S 
CSHT = 6/S 

HT = ARSIN(SNHl) iPM 7094 

I F ( CSHT ) 4 t 5 * 5 

4 HT = FI- H T 
IF(SNHT) 12,5,5 

12 FT - HT-2 ** P l 

5 HAT(J) - FT 

C...HAT - LOCAL TGPOCtNTRIC HUlR ANGLE IN RAD l A NS 
HATGCJ) = H AT ( J ) +XLGNG 

C--..HATG = .GREENWICH TGPG.CENTRIC HOUR ANGLE IN RADIANS 
I POT = IHRUTl J W 
KAT(J) = SDG-HATG ( J ) 

C - . .RAT = TGPOCENTRIC RIGHT ASCENSION IN RADIANS 

WRITE! 6« 214) lOVSIDt I C ) , IDVHN ( J ) , IHRUTC J) , I POT, IRAHG , I RANG , R ASG 
l.RAG, TOCDG, I DC MG , DC SG »DCG , I PRMG , PRSG , PRG * DCT { J > , H A T ( J ) ,HATG( J) , 

2RATU) 

214 FORMAT { 1H , 4 I 3 , 1 X , 1 4 , I 3 , F 7 . 3 , F8 . 5 , I 4 , I 3 , F6 . 2 , F 3. 5 , I 4 , 

1 F 6 . 2 * F 8 * 5 , 3 X , 4 F 8 • 3 ) 

2 CONTINUE 

WRITE 16, 215 ) 

215 FORM AT ( 1H l , 2CX , 4 I 3X , 1 1H rCPUCENTR IC ) / 1 H , 7S : HGC TUBER UT PDT 
1CECUNAT ICN HOUR ANGLE DECLINATION. HOUR ANGLE/ 1H , 74 H DAY 

2 HOUR HOUR DEG PIN DEG MIN DEGREES DEGRE 

3 ES ) 

DO 3 J- 1 , 126 

DC D = DCT( J)*( ISC. /PI ) 

DEC = OCG 
I F ( D CD ) 13,14,14 

13 DEG = -DEC 

14 I one = DEC 
XI DUG = I COG 

XCMN .= ( CEC-aIOUG) *60. 

IF (CCD) 15,16,16 

15 IF (I DUG.) 17,17,18 
17 I DOG = -C 

GO TO 16 



TPCCRu 


cFN 


SOURCE S f A T l i v E N i 


l F »\ t S ) 


10/18/ n 


18 IDOC = - 1 LOG 
16 CONTINUE 

HTO - HATG ( J ) *< 1 80 . /P I ) 
hi = hTD 

IF(HTD) 19,20,20 

19 HT = -HT 
2C 18 ad = hi 

XI H AD a I HA C 

XHPN = ( hi -X I HAL ) *6 C . 

IF(HTD) 21, 2 c t 2 2 

21 IF(IHAL) 23,23,24 

23 I HAD -~0 
GO TO 22 

24 1 H A D = - I Ft AO 

22 CONTINUE 

[POT = IHKUT(J)-7 

toR ITF ( 6, 216 ) I OY MNt J ) , I HKU 1 { J J , I POT , I DOG , XDN-N, IHAC , XHMN, DCD, HTD 

2 16 format.! .1 H , 15 J 1_U , I 6 , I 8 , F 7 . 2 , I 7 , F7, 2 , F 12 .4 , F 14. 41. 

3 CONTINUE 
STOP 
END 
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